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^ ■ ABSTRACT 

^ , Parallel thermal conduction along stochastic magnetic field lines may be reduced 

I because the heat conducting electrons become trapped and detrapped between regions 

fSj ' of strong magnetic field (magnetic mirrors). The problem reduces to a simple but 

, realistic model for diffusion of mono-energetic electrons based on the fact that when 

I there is a reduction of diffusion, it is controlled by a subset of the mirrors, the principle 

^ ■ mirrors. The diffusion reduction can be considered as equivalent to an enhancement of 

. the pitch angle scattering rate. Therefore, in deriving the collision integral, we modify 

f-*) I the pitch angle scattering term. We take into account the full perturbed electron- 

' electron collision integral, as well as the electron-proton collision term. Finally, we 

. obtain the four plasma transport coefficients and the effective thermal conductivity. 

I We express them as reductions from the classical values. We present these reductions 

' as functions of the ratio of the magnetic field decorrelation length to the electron mean 

Q_i' free path at the thermal speed Vr = ^/2kT/me- We briefly discuss an application of 

Q ! our results to clusters of galaxies. 

c/3 . Subject headings: magnetic fields: conduction — magnetic fields: diffusion — methods: 

' analytical — plasmas 
> ' 

^: 

■ 1. Introduction 

The problem of thermal conduction in a stochastic magnetic field is crucial for our under- 
standing of galaxy cluster formation (Suginohara &: Ostriker 1998; Cen & Ostriker 1999) and for 
the theory of cooling flows (Fabian 1990). It is also of great interest for the solar physics and for 
various questions of plasma physics. At the same time, the question: "whether electron thermal 
conduction is so strongly inhibited by a stochastic magnetic field in a galaxy cluster, that it can be 
neglected", is a very controversial one (Rosner & Tucker 1989; Tribble 1989; Tao 1995; Pistinner & 
Shaviv 1996; Chandran &; Cowley 1998). It is currently estimated that if the coefficient of thermal 
conductivity is less than 1/30 of the Spitzer value, then the time scale of the heat conduction in the 
cluster is more than the Hubble time (Suginohara & Ostriker 1998). Otherwise, thermal conduction 
is important. ^ 



^This numerical estimate, 1/30 of the Spitzer value, is based on numerical simulations with limited resolution, so 
it is not the last word on the problem. 



- 2 - 



The problem of thermal diffusion of heat conducting electrons in a stochastic magnetic field 
should be divided into two separate parts because there are two separate effects that reduce dif- 
fusion in the presence of stochastic magnetic field (Pistinner & Shaviv 1996; Chandran, Cowley, 
& Ivanushkina 1999). The first effect is that the heat conducting electrons have to travel along 
tangled magnetic field lines, and as a result, they have to go larger distances between hot and cold 
regions of space. (In other words, the temperature gradients arc weaker along magnetic field lines.) 
The second effect is that electrons, while they are traveling along the field lines, become trapped 
and detrapped between magnetic mirrors (which are regions of strong magnetic field). A trapped 
electron is reflected back and forth between magnetic mirrors until collisions make its pitch angle 
sufficiently small for the electron to escape the magnetic trap. 

In this paper we concentrate on the second effect, and we derive the reduction of the effective 
electron thermal conduction parallel to the magnetic field lines caused by the presence of stochastic 

magnetic mirrors. 

As is well known, a temperature gradient produces electrical current as well as heat flow. 
Similarly, an electric field produces heat flow as well as current. The four transport coefficients 
describing this are given in equation (34) and (35) . The transport coefficients were first calculated 
by Spitzer & Harm for an unmagnetized plasma (Spitzer & Harm 1953; Cohen, Spitzer & Routly 
1950). Their coefficients also apply in an uniform magnetic field for transport parallel to the field. 
In this paper, we show how the parallel transport coefficients can be reduced in the presence of 
stochastic magnetic mirrors, and we calculate their reduced values by the same kinetic approach 
as that of Spitzer & Harm. The reduction factors are presented in Figure 5. The reduced effective 
thermal conductivity (that resulting when the electric field is present to cancel the current) is given 
in Figure 6. Spatial diffusivity of mono-energetic electrons along the magnetic field lines is also 
presented in Figure 3. 

First, in Section 2, we solve the kinetic equation to find the escape time for electrons 
trapped between two equal magnetic mirrors. We assume, that all electrons have a single value 
of speed, V, i.e. they are mono-energetic. The exact calculations of the escape time are given in 
Appendices A and B. In addition, we carry out Monte-Carlo particle simulations to confirm our 
results. 

Second, in Section 3, we apply our results for this escape time to find the reduction of diffusion 
of mono-energetic electrons in a system of stochastic mirrors. It turns out that in the limit Zq S> A, 
where Iq is the magnetic field decorrclation length and A is the electron mean free path, the parallel 
diffusivity is unaffected by magnetic mirrors, and is given by the standard value Dq = {1/3)VX. 
In the opposite limit, Iq <C A, magnetic mirrors do reduce diffusivity. We find that in this case 
there is a subset of the mirrors, the principle mirrors, that inhibits diffusion the most. These 
are mirrors whose separation distances are approximately equal to the electron effective mean free 
path, Aeff , the typical distance that electrons travel in the loss cones before they are scattered out 
of them. In order to estimate the reduction of diffusion in this limit, we need consider only the 
principle mirrors, neglecting all others. Again, we perform the numerical simulations to support 
these theoretical results. 

Third, in Section 4, in order to carry out a precise kinetic treatment involving all electrons, we 
consider the diffusion reduction to be equivalent to an enhancement of the pitch angle scattering 



-3- 



rate of electrons. In deriving the collision integral, we, therefore, modify the pitch angle scattering 
term by the inverse of the factor by which the spatial diffusion is reduced. Wc take into account 
the full perturbed electron-electron collision integral, as well as the electron-proton collision term. 
We obtain an integro-differential equation for the perturbed electron distribution function in the 
presence of stochastic magnetic mirrors. If there is no reduction of electron diffusivity, our equation 
reduces to the well known result obtained by Spitzer and Harm (Spitzer &; Harm 1953; Cohen, 
Spitzer & Routly 1950; Spitzer 1962). 

Fourth, in Section 5, wc solve our equation numerically, separately for the Lorentz gas in 
the presence of magnetic mirrors, neglecting electron-electron collisions (in this case the equation 
simplifies greatly), and for the Spitzer gas in the presence of magnetic mirrors. We find the 
reductions of the four plasma transport coefficients and of the effective thermal conductivity as 
functions of the ratio of the magnetic field dccorrelation length Iq to the electron mean free path 
at the thermal speed Vr = ^/2JcT/me (this mean free path is different for the Lorentz and Spitzer 
models). We find that the major effect of the magnetic mirrors is the reduction of anisotropy of 
superthermal electrons (this anisotropy is driven by a temperature gradient or/and by an electric 
field). Electrical cTirrent and heat are mainly transported by these electrons, whose diffusivity is 
suppressed the most. 

Finally, we discuss our results and give the conclusions in Section 6. 



2. Mono-energetic electrons trapped between two equal magnetic mirrors 

In this section we solve the kinetic equation to find the escape time Tm for electrons trapped 
between two equal magnetic mirrors. We assume here and in the next section that all electrons 
have a single value of speed, V, which is unchanged by collisions, i.e. electrons are mono-energetic. 
In order to derive an analytical solution, we make several additional simplifying assumptions. Let 
the two magnetic barriers (mirrors) be both equal to B,^, and we assume the magnetic field B 

def 

is constant between them. We introduce the mirror strength m = Bt^/B. The separation of 
the mirrors is Z^, and their thicknesses are negligible compared to 1^- In other words, magnetic 
mirrors are similar to thin step-functions with heights — B and with constant field B between 
them (see Figure 1). This is a reasonable assumption, because as we will see in the next section, 
electron diffusion is controlled by strong mirrors with mirror strengths m > 4, which are separated 
by distances much larger than the magnetic field decorrelation length (if the spectrum of mirrors 
falls off with their strength significantly faster than 1/m, the case that we consider in this paper). 

Under these assumptions, the kinetic equation for the distribution function f{t,x,fi) of mono- 
energetic electrons trapped between the two mirrors is (Braginskii 1965) 



df df V d 



2,9 f 



(1) 



Here x is one-dimensional space coordinate along a magnetic flux tube, t is time, ^ = cos 6 is 
the cosine of the electron's pitch angle, and v = V/X is the collision frequency [A is the mean 
free path, see equations (43) and (46)]. The right-hand side of equation (1) represents the pitch 
angle scattering rate, u, of electrons. The electrons are trapped in the region of space between 
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Fig. 1. — (a): A magnetic flux tube with two "step-function like" magnetic mirrors. The mirror 
strengths are m = B^/B. (b): The phase space box where electrons are trapped in coordinates x 
and II = cos 9. The horizontal dotted lines show a closed trajectory of a trapped electron in the 
limit Im *C A. The electrons escape the magnetic trap through two escape windows: x = Imf^, 
M > A*crit = -^1 — 1/m and x = —Im/'^, M < ~/^crit- In the limit Im ^ ^eS the electrons freely 
escape to the right or left whenever they reach the two loss cones, fj, > jUcrit and fj, < — yUcrit- In 
the opposite limit, Aeff ^ Im, the electrons escape when they reach the two shaded regions of the 
phase space. 



the mirrors 
A* > /"crit = \A 



/m/2 < X < Im/'^-, and they can escape through the two windows: 



X 



\/m and x = —Im/'^-, < — A*crit, as shown in Figure 1. The mirror strength is 
"m = B^/B, and it is the measure of the relative heights of the magnetic barriers. For simplicity, 
we assume that the barriers arc high, i.e. m ^ 1 and /Xcrit ~ 1 — l/2m. In this case the electron 
distribution is in quasi-static equilibrium. 



f{t,x,fi) = e ^/^"^F{x,n), Tm:^iy- 



and equation (1) reduces to 



F ^BF V d 
+ fiV— = - — 

Tm ox 2 Ojl 



dF 



(2) 



(3) 



Let us consider an electron traveling in the loss cone /j, > ficrit = ~ 1/m 1 — l/2m (or 

^ < — /Ucrit)- The effective electron mean free path, which is the typical distance the electron travels 
before it is scattered by small angle collisions out of the loss cone, is 

Aeff = A/2m < A. (4) 

In other words, Aeff is a decay distance for a flow of electrons traveling in the loss cones. The 
solution of equation (3) and, therefore, the escape time r^, depends on the mirror strength m and 
the ratio Zm/A. There are three limiting cases for which simple approximate solutions exist: (1) 
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Im <C Aeff = A/2m; (2) Aeff <^ Im A^/Aeff = 2mA; and (3) A^/Aeff <C Im- Wc solve equation (3) 
for case (1) in Appendix A and for cases (2) and (3) in Appendix B, and we obtain the electron 
escape times 

Tm = l^~^lnm, < Acflf, 

T^^ = V-^{lnJKs) =y~^{2mlm/X), Aeflf A V Aeff, (5) 

r(?) = I^-l(3/7^2)(/^/A)^ AVAeff « 

The following simple physical arguments help to understand these results in these three limiting 
cases. The collisional scattering is a two-dimensional random walk of a unit vector (which is the 
direction of the electron velocity) on a surface of a unit-radius sphere with frequency u (so, the 
scattered angle Ag = \/2vt after time interval t). The right hand side of the kinetic equation (1) 
represents a one-dimensional random walk in ^-space that follows from the two-dimensional walk 
because of symmetry. However, it is convenient for the moment to return to the original two- 
dimensional scattering because it is isotropic. The angular sizes of the two loss cones on the 
unit-radius sphere are Aesc ~ First, in the limit Im <^ Aeff, collisions are very weak, and 

the scattered angle over the travel time between mirrors, /m/^, is ~ A/^m/A <C Aesc- Therefore, 
in this case we can disregard the electron motion in x-space. We divide the surface of the unit- 
radius sphere into ~ m boxes, each of angular size ~ Aesc ~ ^/y/rn. The time it takes for the 
unit vector to random walk from one box to another is ~ i'~^/m, resulting in the total escape 
time Tm fn X {u~^ /m) = u^^. Because the unit vector can "visit" each box more than once, 
the exact result contains the logarithm of m. Second, in the limit Aeff <^ Im ^ A^/Aeg, wc have 
to consider motion in x-space as well. In this case the electrons move in three-dimensional phase 
space, and they escape when they are in the two loss cones within distance Aeff from the mirrors, 
as shown by the shaded regions in Figure 1(b). We divide the three-dimensional phase space into 
~ (Zm/Acff)(l/A^gj,) ~ m^lm/X boxes, each of size Aeff A^g^, ~ X/m? . The time it takes to move from 
one box to another is ~ u~^/m, resulting in the total escape time ~ (m^i^/A) x (z/~^/m) = 
v~^{mlm/X). Note, that the electron distribution function is almost constant in the phase space in 
this case (see Appendix B). Third, in the limit A^/Aeff <C Im-, the escape of electrons is controlled 
by slow diffusion in x-space, so the escape time is approximately equal to the time of diffusion 
between mirrors, ~ '^^^(^m/A)^ in this case. 

In our further calculations we use a simple interpolation formula 



^(1) + ^(2) + .^(3) = J,- 
'm ^ 'm ^ 'm 



Inm + {lm/\w) + (3/7r2)(Z^/A)' (6) 



for the whole range of parameters m and Im/X. This formula is suggested by the numerical sim- 
ulations shown in Figure 2. The dots in this figure show the results of our Monte-Carlo particle 
simulations for three mirror strengths m = 2, m = 16 and m = 128. To obtain these results we 
followed 10^-10^ electrons trapped between two equal magnetic mirrors separated by distance 
ranging from 1/1024 to 256 in units of the mean free path A. Independently of the initial distri- 
bution of electrons, the number of trapped electrons tends to an exponential dependence on time 
with the characteristic decay time Tm in just a few collision times [see equation (2)]. The solid lines 
in the figure represent formula (6) and are in a very good agreement with the simulations even for 
the smallest mirror strength m = 2. 
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Fig. 2. — The dots show a logarithmic plot of the numerically obtained electron escape time 
in units of the collision time as a function of the separation 1^ of two equal magnetic mirrors 

in units of the mean free path A. These results are based on our Monte-Carlo particle simulations 
of 10^-10^ trapped electrons, assuming three values of the mirror strengths, m = 2, m = 16 and 
m = 128. The solid lines represent the analytical result, equation (6). 

3. Diffusion of mono-energetic electrons in a system of random magnetic mirrors 

In this section we continue to assume that electrons have a single value of speed, V. If there were 
no magnetic mirrors and the magnetic field had constant strength along the field lines, the parallel 
diffusion of mono-energetic electrons would be the standard spatial diffusion, Dq = (1/3)FA. Here, 
A is the electron mean free path at speed V. However, as wc have discussed in the introductory 
section, diffusing electrons move along flux tubes with random magnetic field strength and become 
trapped and detrapped between magnetic mirrors. These mirrors are regions of strong field and 
are separated by a field decorrelation length Iq. As a result, the diffusion is reduced by a factor 
that depends on the ratio lo/X. 

In the main part of this section we derive this diffusion analytically and at the end of the 

section confirm it with numerical simulations. (In contrast to the previous section, where there 
were only two equal mirrors, in this section, we consider many mirrors with random spacing and 
strength.) 

Consider the limit ^ first. In this case collisions are strong, and according to the third 
formula in equation (5), the time it takes for electrons to escape a trap between two magnetic 
mirrors is independent of the mirror strengths and is entirely controlled by the standard spatial 
diffusion transport of electrons between the mirrors. As a result, magnetic mirrors can be ignored, 
and there is no reduction of diffusion; D = Dq. 

In the opposite limit, Iq <^ A, the collisions are weak, and magnetic mirrors do result in a 
reduction of diffusion. To find this reduction, we divide all mirrors into equal size bins bm = 
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(m — 5/2, m + (5/2], where m is the bin central mirror strength, and constant 6 is the width of the 
bins (the value of 6 will be discussed later) . 

For the moment we consider the diffusion in the presence of only those mirrors that are in a 
single bin 6^. It turns out that one of the bins leads to a smaller diffusion than any other bin, 
and the net diffusion due to all the mirrors is approximately that due to only mirrors in this bin, 
provided that the bins are sufficiently wide. 

Let the spectrum of magnetic mirror strengths be V{m). We assume that strong magnetic 
mirrors are rare, i.e. the spectrum falls off fast with the mirror strength (we will estimate how fast 
it should fall off, below). The probability that a mirror belongs to bin hm is 

Pm = I V{m')dm' K5V{m) + {5"^/2A)V"{m). (7) 

At each decorrelation length Iq the magnetic field changes and becomes decorrelated. Therefore, 
the mean separation of mirrors that are in bin hm is 

oo 

lm = h^kpra{l-Pmf~^ = hlPm- (8) 
fe=l 

Let us consider an electron trapped between two mirrors of bin hm- The time Tm that it takes 
for this electron to escape the trap is given by equation (6), where we keep only the first two terms 
(because Iq <^ A) 

Tm ~ t!^^ + =iy~Hn {mqm) . (9) 
Here, we introduce the important parameter 

Qm =^ exp (Im/Ks) = exp (2mlo/PmA), (10) 

where the mean distance between the two mirrors is given by equation (8). After the electron 
escapes, it travels freely in the loss cone in one of the two directions along the magnetic field 
lines until it is again trapped between another two mirrors of bin 6^. The freely traveling electron 
becomes first trapped with probabilities 1 — e~''"/'*'«ff = 1 — in < x < Z^, e~'"'/'^«ff — e"^'"'/'^"'*^ = 
1m^~1m^ in < X < 2lrn, e"2''"/'^<=ff - e~^''"/'^<=ff = q~i^-qm^ in 2l„i < x < 3lm, and so on. Therefore, 
the mean distance squared (Ax^)^ that the electron travels in the loss cones before trapping is 

(Ax^)^ - 4 f; fe^(g^'=+^ - .^^i = 4 7^^^^- (11) 
k=i ^^"^ ^> 

The processes of trapping and detrapping repeat in time intervals r„j. In other words, electrons 
random walk along the field lines in a system of mirrors that belong to bin bm with steps ~ 
(Ax^)^ in time intervals « Tm- As a result, the diffusion coefficient for these electrons is D{m) = 
C [(Ax^)^/2rTO], where we introduce a scaling constant C, which is of the order unity and will be 
determined by the numerical simulations. The corresponding reduction of diffusion is 

n/ ^/n ^3 ( gm(gm + 1) 1 1 , 

0(m)/D„ = C-(^-j i.«A. (12) 
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where we use Do = (l/3)z^A^ and equations (8), (9) and (11); Pm and qm are given by equations (7) 
and (10). 

For a given spectrum of mirrors V{m) and given constants 5 and C, the diffusion reduction (12) 
due to mirrors of bin b^n, is a function of mirror strength m. Let us analyze this function in two 
hmits: Ing^, *C 1 and lng„ 3> Inm > 1. If Ing^ <C 1, then — 1 = 2mlo/pmX <C 1. Therefore, 
D{m)/Do C(3/4)(l/m2 Inm) and {d/dm)[D{m)/DQ] < 0. On the other hand, if Ing^ » Inm, 
then D{m)/Do C(3/4)(;o/A)(l/m;),„). Therefore, (d/(im)[L>(m)/I?o] > if the spectrum of 
mirrors falls off faster than 1/m with the mirror strength.^ In this paper we make an assumption 
that the spectrum falls off significantly faster than 1/m. 

Therefore, a minimum of D{m)/Do exists. Let this minimum be achieved at m = nip. Then 
^^Qrup = Irup/^eS ~ 2/lnmp ~ 1, or Imp ~ Aeff. The minimum can roughly be estimated as 
D{mp)/DQ = mm {D{m)/ Do} ~ ^/rrip, which is in agreement with the qualitative results of Al- 
bright et al. (2000). 

In other words, if Iq <^ A, then there is the bin that inhibits diffusion the most. We call it the 
principle bin, bp = {nip — 6/2, nip + S/2]. The corresponding mirror strength nip is the principle 
mirror strength. The minimum of diffusion D{ni) due to mirrors of bin bm is achieved at the 
principle strength, m = nip. The spacing of mirrors that are in the principle bin is of the order of 
the effective mean free path for this bin, ~ Aeff = \/2nip. The main idea is that, in order to 
estimate the net diffusion due to all mirrors, we need consider only magnetic mirrors that are in the 
principle bin and we can neglect all other bins. Mirrors that are smaller than the principle mirrors 
"work" poorly in the inhibition of diffusion because they are weak and are separated by distances 
less than Aeff (which is the distance that electrons travel in the loss cones). Mirrors that are larger 
than the principle mirrors "work" poorly, because they are very rare and are separated by very 
large distances (provided the mirror spectrum falls off with the mirror strength significantly faster 
than 1/m). These assumptions are supported by our numerical simulations (see Figure 3). 

As a result of these considerations, we can combine our theoretical results for the reduction of 
diffusion of mono-energetic electrons, i?D = D /Do, into a single formula valid in the two limits for 



where D{ni)/Do is given by equation (12), and the minimum is achieved at the principle mirror 
strength m = nip (note that Ing^p = Inip/^cS ~ 1)- 

We show the theoretical mono-energetic diffusion reduction (13) by the solid lines in Figure 3 



^This criterion is different from the result of Albright et al, (2000), who found to be the boundary spectrum 

for the transition between their diffusive and subdiffusion regimes. We believe that the difference arises because, for 
flat spectra, our bin width S starts to depend on lo/X (and our simple diffusion model breaks down). 



Zo/A: 




h < A, 
lo » A, 



(13) 
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(a) Exponential: C=2.13, (5=2.60 (b) Gaussian: C=1.68, (5=1.07 




Fig. 3. — We consider two mirror spectra: (a) exponential, and (b) Gaussian [see eqs. (14)]. The 
dots show the reduction of diffusion, = D/Dq, obtained by Monte-Carlo particle simulations 
of 1-6 X 10^ electrons, each followed in a system of magnetic mirrors over 300 collision times u~^. 
The solid lines represent the theoretical results given by equation (13). The constants C and S 
are obtained by matching the theoretical results with the results of simulations for each of the two 
spectra (and these constants do not depend on Iq/X). 



(a) Exponential: Zo/A=l/l6 



(b) Gaussian: lQ/\=^/^6 
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Fig. 4. — The natural logarithm of the diffusion reduction (12) caused by mirrors that are in bin bm 
for /o/A = 1/16. We consider two mirror spectra: (a) exponential, and (b) Gaussian [see eqs. (14)]. 
The principle bins are shown by arrows. In case of each spectrum, the reduction has the minimum 
at the principle mirror strength nip, and it roughly doubles at the boundaries of the principle bin. 
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for two mirror spectra: exponential and Gaussian^, 



P(m) = e ^'^ exponential, 

V{m) = (2/7r)^/^e-("*-2)V2 Gaussian. 



(14) 



The results of our Monte-Carlo particle simulations are shown by dots. The constants C and 6 
(shown at the top) are of the order unity, and we adjust them by matching our theoretical results 
with the results of simulations in case of each of the two spectra (C and S do not depend on Iq/X). 
The simulations are based on 1-6 x 10^ particles. For each particle we choose a distribution of 
mirrors m > 2, which all are separated by the magnetic field decorrelation length and are chosen 
according to the assumed mirror spectrum (14). We follow the particles during 300 collision times 
u~^. Then wc average the particle displacements squared (Ax^) at a given time t to obtain the 
diffusion coefficient {Ax'^)/2t given in Figure 3. 

Note that the bin width 5 is larger for the exponential spectrum than it is for the Gaussian. This 
is because the later is steeper at large mirror strengths. Figures 4(a) and 4(b) clearly demonstrate 
the difference. In these figures we plot the natural logarithm of the diffusion reduction (12) caused 
by mirrors that are in bin bm versus the mirror strength m for Iq/X = 1/16 and for the both 
spectra (14) of mirror strengths. The principle bins are shown by arrows. In case of each spectrum, 
the reduction has the minimum at the corresponding principle mirror strength nip. We sec, that 
the reduction roughly doubles over its minimal value at the boundaries of the principle bin, m = 
nip + 5/2 and m = nip — 6/2. 



In this section we use the results found above to obtain the modified kinetic equation for 
electrons traveling in a system of random magnetic mirrors. The reduction of diffusion of mono- 
energetic electrons with speed V, Rn, obtained in the previous section can be considered to be 
equivalent to an enhancement of the pitch angle scattering rate, since the pitch angle scattering 
is directly related to spatial diffusion. We therefore, in deriving the collision integral, modify the 
pitch angle scattering term by factor R^^ , where Rd is the factor by which the spatial diffusion is 
reduced (see the previous section). Hereafter, we do not assume electrons to be mono-energetic. 
We take into account the full perturbed electron-electron collision integral, as well as the electron- 
proton collision term. When i?D = Ij our equations reduce to those of Spitzer and Harm in their 
well known paper (Spitzer & Harm 1953; Cohen, Spitzer & Routly 1950; Spitzer 1962). 

The electron distribution function is 



4. The Fokker-Planck kinetic equation 



f{^l,v) = foiv) + h{^i,v) 



(15) 



where /o is the zero order isotropic part given by the Maxwellian distribution 



/o = n{x) [nie/27rkT{x)f^ 




/2kT{x) ^ 



(16) 



^We find the minimum in equation (13) numerically. 
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and fi (X fxis the first order anisotropic perturbation (of order the temperature gradient and electric 
field) 

fiifi,V) = i^nV^^S{v). (17) 

Here nie is the electron mass, k is the Boltzmann constant, and the electron temperature T(x) 
and concentration n{x) slowly change in space. We also introduce the dimensionless electron speed 
V = V/Vt, where the thermal electron speed is Vt = \J1kT jm^,. Thus, the function S'(f) in 
equation (17) is dimensionless. 

In a steady state, the kinetic equation for the electrons is obviously 

V^dfo/dx) - {eE/me){dh/dV:,) = {Sf/St)^, (18) 

where {Sf /St)^ is the Coulomb collision integral that includes cicctron-proton and electron-electron 
collisions, Vx = fiV is the x-component of the electron velocity (the component along the magnetic 
field lines), and E is the electric field in the x-direction. The electron pressure should be constant, 
P = kn{x)T(x) = const. ^ As a result, the derivatives of the Maxwellian electron distribution are 

dfo/dx = {v''-2.5){fo/T){dT/dx), dfo/dVx = -{2n/VT)vfo. (19) 



The collision integral is divided up as 



(20) 



where (5/o/(5t)g = corresponds to Maxwellian collisions acting on /o, {Sfi/St)^ corresponds to 
Maxwellian collisions (with enhanced pitch angle scattering) acting on /i, and {Sfo/6t)^ corresponds 
to perturbed collisions acting on /o (since /o is isotropic, there is no pitch angle scattering in 
this collision term). The collision integral (20) can be best obtained, in the Fokker-Planck form, 
by using the Rosenbluth potentials h{p,V) = hQ{V) + hi{fi,V) and g(iJ,,V) = go{V) + gi{p,V) 
(Rosenbluth, MacDonald, & Judd 1957). Here ho and go are calculated using the Maxwellian parts 
of the electron and ion distribution functions (16), while the perturbed potentials, hi = 2iJ,Ai{V) 
and gi = fiBi{V), are proportional to /x, and they are calculated using the perturbed part of the 
electron distribution function (17). 

The Maxwellian potentials ho and go determine the {dfi/dt)^ part of the Fokker-Planck col- 
lision integral, and the perturbed potentials, h\ = 2nAi{V) and g\ = jiBiiV), are used to find 
the ((5/o/(5t)i part of the Fokker-Planck collision integral [see the equation (31) of Rosenbluth, 
MacDonald, & Judd 1957] 



Ad 
2n 

J__d_ 
y2 gy 



LA 
V^dV 

dgo 



2 dho 



h 



dV 



•'^ dV 



+ 



1 

2y2 

dgo d 



2V^ dV df^ 



fiV 



dv^ 

dfi 



(21) 



''Because the hydrodynamic time scale is much shorter than the transport, e.g thermal conduction, time scale. 
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Ad 
2n 



2 d 
V^dV 



foV' 



,dAi 
'dV 



foV 



For a hydrogen plasma the "diffusion constant" Ad is 



_3__ci_ 
"^V^dV 



Ad = Sirne^ In A/ ml 



(22) 



(23) 



where e is the absolute value of the electron charge, and In A is the Coulomb logarithm (Spitzer 
1962). Note, that the last term in equation (21) is the pitch angle scattering term, and we multiply 
it by our enhancement factor R^^ [compare this term with the right-hand side of equation (1)]. 

Using the equations (17) and (18) of Roscnbluth, MacDonald, & Judd (1957), we express 
the derivatives of the potentials ho and go in terms of the three Maxwellian diffusion coefficients 
(AV||)q, {{AV±)'^)q and ((AV||)^)q, which are further given in terms of error functions [see the 
equations (5-15)-(5-20) of Spitzer 1962] 



dho/dV 
dgo/dV 



{2n/An) (AV||)q = -{n/V^)[l + 4v^G{v)], 



{n/AD)V ((AVj_ 



n[l + ^{v) - G{v)], 



(24) 



{2n/AD){iA^r)^ = {2n/V)G{v). 



Here $ is the usual error function, and G is expressed in terms of $ and its derivative they are 
functions of the dimensionless speed v = V/Vt [Vt = y^fcT/ml], 



= (2/v^) 

Jo 



e ^ dx, G{v) = 



2^ ■ 



(25) 



The perturbed potentials, hi = 2fiAi{V) and gi = fjBi{V), arc calculated using the perturbed 
electron distribution function (17) and are given by the following formulas [see the equations (40), 
(41), (45) and (46) of Rosenbluth, MacDonald, & Judd 1957] 



^1 = iA7r/3){n/VT)[v-^l3{S;v)+vlQiS;v)], 

Bi = (47r/3)n1/T [^.2v-^7^{S; v) - lz{S- v) - vl2{S; v) + 0.2v^lo{S; v) 



where we introduce integrals 

'lm{S;v) 



V^S{v)dv, Lm{S;v) 



f 



v'^Siv) dv, 



(26) 



(27) 



Now, substituting equations (16), (17), (24) and (26) into formulas (21) and (22), and using 
definitions (25), (27) and equation (20), after considerable algebra, we have for the collision integrals 

(<5/i/5t)o = {nAD/2VS)iJv-\CS-2v^<!>'S), 
{6fo/6t)-^ = {nAD/2V^)fiv~^{iS + 2v^<^'S), 

(Sf/St)^ = {nAD/2VS)fiv-^iCS + IS), (28) 
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where the differential and the integral operators are defined as 

jC,S{v) = d/dv[vG{dS/dv)]+2v'^G{dS/dv)-[v-^R^\l + ^-G)-4v^^']S, (29) 
iS{v) = (4/150F) e~"' [12/5(5; v) - 10/3(5; v) + 2v^{6v'^ - 5)Iq{S; v)] . (30) 

The enhancement of the Maxwellian pitch angle scattering rate, Rji)^ , enters into the differential 
operator (29). Rb depends on the dimensionless speed v = V/Vr, we will explicitly give this 
dependence in equations (42) and (45). 

Finally, substituting formulas (28) and (19) into equation (18), we obtain the kinetic equation 
for the dimensionless perturbed electron distribution function S{v) [see equation (17)] 

Is = 7t v^{2v^ - 5)6"'^' + 7e v^e-""^ - IS, (31) 
S{v) — 0, as f — and as f — 00, (32) 

where constants 7t and 7b are 

k^T dT _ kT 

27rV2ne4lnA d^' ~ ix^l^ne^hiK ' ^ ' 

Wc also take the obvious boundary conditions (32) for function S. Equations (29)-(31) reduce to 
the Spitzer equations for an ionized hydrogen gas (Spitzer &; Harm 1953; Cohen, Spitzer k. Routly 
1950) if we set R^, = 1 and make a substitution S{v) = 7r'~^/^e~'^ D{v). However, we prefer to use 
function S, because of the simpler boundary conditions (32). 



5. The reduction of transport coefficients by stochastic magnetic mirrors 

In a steady state, an electric field E and a temperature gradient dT/dx both produce anisotro- 
pic perturbations of the electron distribution function, fi{fi,v) = finVrj^^ S (v) , see equations (15) 
and (17). This anisotropy results in an electron flow and, consequently, in an electric current j and 
in a heat flow Q along magnetic field lines (in the x-direction) 

j = -e / iiV fidn2TrV^dV = aE + a (dT/dx), (34) 

poo pi 

Q = J J nV{meV^/2)fidn27rV^dV = -PE-K (dT/dx). (35) 

Here a, a, (5 and n are the four transport coefficients to be found (a and k are the electrical and 
thermal conductivities). 

Before we proceed to the calculation of the transport coefficients, let us first call attention to 
the electron flow produced by the electric field. The electric field produces two different kinds of 
the electron flow. The flrst, the main, flow is due to acceleration of electrons, which is described 
by the term containing E in equation (18), and correspondingly by the term containing 7e in 
equation (31). The second, an additional, flow arises because the electric field changes the size of 
the two loss cones of a mirror trap, so in Figure 1(b) /Xcrit in the right upper corner is not equal to 
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Hcrit in the left lower corner. As a result, the electrons are more likely to escape the trap in the 
direction opposite to the electric field. Fortunately, this additional flow, which is rather complicated 
to find precisely, can be neglected compared to the fiow due to acceleration. We give a prove of 
this in Appendix C. ^ 

In further calculations, it is convenient to break S(v) into the two separate inhomogeneous 
solutions of equation (31), which we denote as St{v) and Se{v). ^ The first solution, St, is 
obtained by setting 7t = 1 and 7e = 0, and the second solution, Se, is obtained by setting 7t = 
and 7b = 1, i. e. 

St{v) = S{ij), when 7x = 1 and 7e = 0, 
SEiij) = S{v), when 7t = and 7e = 1. 

The general solution to equation (31) and the perturbed distribution function (17) are the linear 
combinations of the two inhomogeneous solutions, 

S{v) = -fTSTiv)+-fESE{v), 
fi{^i,v) = iinV^^[j^STiv)+j^SE{v)]. (37) 

In other words, St and Se correspond to anisotropic perturbations of the electron distribution 
function, which are driven by the temperature gradient and by the electric field respectively, while 
S = 7tS't + JeSe is the total anisotropic perturbation. 

We now consider separately two cases: first, the Lorcntz gas in a system of random mirrors, 
and second, the Spitzer gas in a system of random mirrors. For the Lorentz gas, electrons are 
assumed only to collide with protons, so equations (29)-(31) become greatly simplified. For the 
Spitzer gas, we consider both the electron-electron the electron-proton collisions, so we solve the 
full set of our equations. 



5.1. Lorentz gas in a system of random mirrors 

Here we assume the electrons to collide only with protons, so we have for operators (29) 
and (30) 

CS = -S/vRb, is = 0, (38) 
resulting in the two simple inhomogeneous solutions (36) of equation (31), 

St{v) = -v\2v'^ - 5) e-"'i?D, Se{v) = -v^ e""' Rb- (39) 



^The main reason is that the difference in the two loss cones due to electric field is inversely proportional to the 
electron kinetic energy, so the additional flow has a factor 1 /V'^ compared to a factor 1 /V^ that enters the main flow 

due to acceleration. Because both the current and the heat flow are mainly transported by superthermal electrons 
V — V/Vt ~ 2, the additional flow is approximately 20% of the main flow, see Appendix C. 

®The two homogeneous solutions of equation (31) must be excluded, because they diverge either at u — » or at 
w — > oo, violating the conditions (32), see more details in Cohen, Spitzer & Routly 1950. 



(36) 
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If there are no magnetic mirrors, so i?D = 1, wc substitute equations (39) into formula (37) and eas- 
ily carry out the two integrals in equations (34) and (35). Taking into consideration definitions (33), 
we obtain the well-known Lorentz transport coefficients (Spitzer 1962) 

o-L = 2 - wo _ > aL = 3 - 



''^ / me^^e^lnA x''^ ) me'^^e^lnA 

— 1 ttt: , Kt — 20 I — 1 T-r: 



'^^'[-^) ^^^'%) 



If there are magnetic mirrors, it is convenient to normalize the resulting transport coefficients 
to the corresponding Lorentz coefficients (40). Substituting equation (37) into the two integrals in 
equations (34) and (35), and again using definitions (33), we have 

cr/(7L = -(1/3)T3(5e;oo), a/a^^ = -(1/9) T3(5t; oo), 

;3//3l = -(1/12)75(5e;oo), «K = -(1/60)75 (5t;oo), (41) 

where the integral moments are defined by equations (27), and St and Se are given by equa- 
tions (39). 

In order to find explicitly the diffusion reduction factor Rd in equations (39) as a function 
of u, we refer to the results of Section 3. In those section we found the diffusion reduction as a 
function of the ratio of the magnetic field decorrelation length Iq to the electron mean free path 
A. For Lorentz electrons the mean free path Al is proportional to the fourth power of the electron 
speed, Al oc F^, (Spitzer 1962, Braginskii 1965). Thus, we have 

i?D = Rb{Io/Xl) = Rb{v-%/Xl,t), (42) 
where Al,t is obviously the Lorentz electron mean free path at the thermal speed Vr = y^2kT/me 

Al,t= (A;r)7vrne^lnA« 0.1Kpc(T/10^K)^(10-2cm-Vn). (43) 

Here we assume the Coulomb logarithm for a cluster of galaxy to be InA 40 (Suginohara & 
Ostriker 1998). 

Wc use our theoretical results given by equation (13) for the mono-energetic diffusion reduction 
Rd = Rd{1'o/Xl) = Rd{v^'^Io/ Xl,t) in the limits w^^/o/Al,t ^ 1 and v~'^Iq/ X-^^t ^ 1; and we use 
our numerical simulation results presented in Figure 3 for f~'^/o/-^L,r ~ 1- [We carry out the cubic 
spline interpolation of the simulation results. Note, that i?D is not differentiated in operator (29), 
so our final results are not sensitive to small noise errors in calculation of i?D-] 

Using equations (39) and (42) with i?D given in Section 3, and numerically performing the 
velocity integrals, we find all four transport coefficients (41) normalized to the standard Lorentz 
coefficients (40). The dashed lines in Figures 5(a)-(h) show the resulting normalized transport 
coefficients a, a, (5 and k as functions of Io/X^^t for the two mirror spectra: (a) exponential, and 
(b) Gaussian [see equations (14)]. The asymptotic values of the coefficients at large values of Zo/Al,t 
are given by the numbers on the dashed lines, and they are unity. Thus, there is no reductions 
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Fig. 5. — Figures on the left/right correspond to the exponential/Gaussian mirror spectra [see 
eqs. (14)]. The sohd/dashed hnes on all figures show the four transport coefficients, a, a, (3 
and K, for the Spitzer/Lorcntz gas in the presence of stochastic magnetic mirrors as functions of 
the ratio of the magnetic field dccorrelation length /q to the Spitzer/Lorcntz electron mean free 
path, As,t/Al,t, calculated at the electron thermal speed Vt = y2fcr/m^ [see eqs. (43), (46)]. 
All transport coefficients are normalized to the standard Lorentz transport coefficients given by 
equations (40). The asymptotic values of the coefficients at Iq/Xt 3> 1 are given by the numbers 
on the lines. They agree with the results of Spitzer &; Harm (1953). 
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of the transport coefficients at Io/^l,t ^ 1, as one can expect because tfiere is no reduction of 
electron diffusivity in this hmit [see equation (13)]. 

In a steady state, the electrical current j in a highly ionized plasma should be zero. Thus, if 
a temperature gradient is present, the resulting electric field E is obtained by setting j to zero in 
equation (34). Substituting this result for E into equation for the heat flow (35), we find for the 
effective thermal conductivity 

Kes = K — af3 1 a, 
Kcff//«L = n/n^- i3/5)ia/a^){P/P^){aJa), (44) 

where we use formulas (40) for the Lorentz transport coefficients in the second line of this equation. 

Using the transport coefficients reported in Figures 5 by dashed lines and formula (44), it is 
easy to find the effective thermal conductivity k^s normalized to the standard Lorentz thermal 
conductivity [see equation (40)]. However, it is more useful to give the ratio of k^s to the 
Lorentz effective conductivity, /«L,efr = 0.4 This ratio is the actual suppression of the effective 
conductivity of the Lorentz gas by magnetic mirrors. The dashed lines in Figures 6 show this 
suppression, Kcflf/^L.cff) functions of Io/\ij,t for the two mirror spectra: (a) exponential, and (b) 
Gaussian [see equations (14)]. It has been estimated that the time of heat conduction in clusters of 
galaxies is possibly larger than the Hubble time if the thermal conductivity is less than 1/30 of the 
Spitzer value (Suginohara &; Ostriker 1998). The horizontal dotted lines indicate this reduction of 
1/30. 

For comparison, the dotted lines represent the mono-energetic diffusion reduction at the elec- 
tron thermal speed, Rd{Io/^l,t) = D{Iq/ Xi^^t)/ Dq. We see, that the Lorentz gas effective conduc- 
tivity is reduced to a value two to three times smaller than that of the diffusion reduction. This 
is because heat is mainly transported by superthermal electrons. These electrons have long mean 
free paths, and the magnetic mirrors more strongly inhibit their diffusion. 



5.2. Spitzer gas in a system of random mirrors 

Now consider the full collision integral (28) for the Spitzer gas in a system of random magnetic 
mirrors. We have numerically solved the full set of our equations (29)-(32). The formulas (41) 
and (44) remain the same as for the Lorentz gas, but the functions St{v) and Se{v) are different. 
For Spitzer electrons the mean free path is As oc y^[l + ^(v) — G{v)]~^ [Spitzer 1962, the error 
functions $ and G are given by (25)]. Thus, formula (42) for the reduction of spatial diffusivity 
now becomes 

where the Spitzer electron mean free at the thermal speed Vt = ^JlkTjme is 

As,T = 0.614 (A;r)77rne^lnAf« 0.06Kpc(r/10^K)^(10"3cm-V")- (46) 
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Fig. 6. — The solid/dashed lines show the reduction of the parallel effective thermal conductivity 
for the Spitzer/Lorentz gas by stochastic magnetic mirrors, as function of the ratio of the magnetic 
field decorrelation length to the Spitzer/Lorentz electron mean free path. The notations are the 
same as in Figure 5. For comparison, we give the mono-energetic diffusion reduction Rj:){Iq/ Xt) = 
D{Iq/ Xt)/ Dq by the dotted lines. The horizontal dotted lines represent the reduction of 1/30, 
below these lines the thermal conduction is so weak, that it should become negligible in clusters of 
galaxies. 
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Fig. 7. — The solid lines show functions v^St{v) [the left figure] and v^Se{v) [the right figure] for 
the Spitzcr gas in a system of random magnetic mirrors for the case Iq = As,t [^o is the magnetic 
field decorrelation length, As,r is the Spitzer electron mean free path (46)]. The dashed lines in the 
corresponding plots show the same functions for the Spitzer gas without mirrors. Both graphs are 
plotted for the exponential mirror spectrum (for the Gaussian spectrum the results are similar) . 
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Functions St{v) and SEiv) arc defined by equations (36), and they are the two inhomogeneous 
solutions of the equations (31), (32). To find these solutions we solved equation (31) numerically 
by iterations. At each iteration step the integral part of this equation, IS, was calculated using the 
solution for S from the previous step, and the new solution for S was calculated by the Gaussian 
decomposition with backsubstitution (Fcdorcnko 1994), using the boundary conditions (32). Ini- 
tially, wc started with zero function = 0. The iterations converged very rapidly, and the Gaussian 
decomposition method is stable. We believe that our numerical method is much better and faster 
than the method of Spitzer and Harm (1953) because their method was not stable. It took us less 
than ten seconds of computer time to calculate all digits of the transport coefficients reported by 
Spitzer and Harm. 

The solid lines in Figures 5(a)-(h) show the resulting transport coefficients a. a, (3 and k 
normalized to the standard Lorentz coefficients (40) as functions of Iq/Xs.t for the two mirror 
spectra: (a) exponential, and (b) Gaussian [see equations (14); remember that /q is the magnetic 
field decorrelation length] . The asymptotic values of the coefficients at large values of Iq / As,r are 
given by the numbers on the solid lines, and they agree with the results of Spitzer and Harm. 

The effective thermal conductivity, Keff, normalized to the Spitzer effective conductivity, Ks,eff = 
0.0943 Kl, is given in Figures 6 by the solid lines for the two mirror spectra. This normalized 
conductivity is the actual suppression of the effective thermal conductivity of the Spitzer gas by 
stochastic magnetic mirrors. It is the result that should be applied in astrophysical problems with 
random magnetic mirrors. 

Finally, it is interesting to see how the mirrors change the Spitzer perturbed electron distri- 
bution function. In Figure 7 we plot functions v^St{v) and v'^Se{v) for the case when Iq = As,t 
[note that 2ttV'^ f\ is the actual distribution of electrons over speed V = vVr, see equation (37)]. 
The solid lines represent these functions for the Spitzer gas in a system of random mirrors with 
the exponential mirror spectrum (for the Gaussian spectrum the results are similar). The dashed 
lines show the same functions for the Spitzer gas without magnetic mirrors. Wc sec that v^St{v) 
and v'^Se{v) are reduced at large values of i.e. magnetic mirrors reduce the anisotropy of the 
superthermal electrons, which carry the electrical current and heat. 



6. Conclusions 

In this paper we have derived the actual parallel effective thermal conductivity that should 
be applied to astrophysical systems with random magnetic mirrors, as well as other important 
transport coefficients. 

Now, let us apply our results for the reduction of the Spitzer effective electron thermal con- 
ductivity, shown in Figure 6 by the solid lines, to the galaxy cluster formation problem. If the 
reduction is by more than a factor of thirty (shown by the horizontal dotted lines in Figure 6) , then 
the time of heat transport becomes larger than the Hubble time, and the heat conduction can be 
neglected (Suginohara k, Ostriker 1998). ^ We see that this is the case if the magnetic field decorre- 



See the footnote on page 1. 
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lation length Iq is roughly less than 10^^ - 10^^ of the electron mean free path at the thermal speed 
At = \{\/2kT /rrig) (we consider the Spitzer gas). Although there is little observational data about 
the topology of magnetic fields in clusters of galaxy, the magnetic field scale is probably 1-10 Kpc 
(Kronberg 1994; Eilek 1999 and references in it). According to equation (46) the characteristic 
electron mean free path at the thermal speed is 0.06 -60 Kpc for temperatures T = 10^-10* K and 
densities n = 10^^-10"^ cm~^. We see, that in general, the effective electron thermal conductivity 
parallel to the magnetic field lines is not reduced enough by magnetic mirrors to be completely 
neglected. However, as we pointed out in the introductory section, there is an additional effect 
that electrons have to travel along tangled magnetic field lines larger distances from hot to cold 
regions of space, so the thermal conduction is further reduced (this effect will be considered in our 
future paper). At the moment, "whether electron thermal conductivity in clusters of galaxies is 
sufficiently inhibited that it can be ignored" is still an open question. 

Recently, Cowley, Chandran et al. studied the reduction of the parallel thermal conduction, 
and they concluded that the thermal conductivity in galaxy clusters is reduced enough to be 
neglected (Chandran k. Cowley 1998; Chandran, Cowley k, Ivanushkina 1999; Albright et al. 2000). 
Their conclusions are different from ours. The reason is that our approach in calculation of the 
conductivity is very different, and our results are qualitatively different. The main difference is 
that they took the reduction of thermal conductivity to be equal to the reduction of diffusivity of 
thermal electrons. In fact, the reduction of diffusivity is due to the enhanced pitch angle scattering 
by stochastic magnetic mirrors, and to find the reduction of thermal conductivity, the full set of 
kinetic equations must be derived and solved. This consistent way of solving the problem makes a 
considerable difference (see Figure 6). On the other hand, Cowley, Chandran and et al. first called 
attention to the importance of the effective mean free path Aeff and found the correct qualitative 
result, that in the limit Iq ^ X the diffusion reduction is controlled by the mirrors whose spacing 
is of order of the effective mean free. 

We are happy to acknowledge many useful discussions of this problem with Jeremiah Ostriker, 
Jeremy Goodman and David Spergel. We would also like to thank Makoto Matsumoto, Takuji 
Nishimura and Shawn J. Cokus providing us with fast random number generators (which are given 
at http:/ /www. math. keio.ac.jp/matumoto/emt.html). This work was partially supported by DOE 
under Contract No. DE-AC 02-76-CHO-3073. Leonid Malyshkin would also like to thank the 
Department of Astrophysical Sciences at Princeton University for financial support. 



A. Solution of equation (3) in the limit Im *C Aeff 

Here we solve equation (3) by expansion in the limit ^ Acfj. This condition means that 
collisions are too weak to scatter the electron out of the loss cones. Therefore, F(x,ii) = when 

We make use the fact that {V /u){d / dx) ~ X/lm 3> 1. Also we will show that l/vTm <C 1. 
The validity of this last assumption appears below. To zero order, we have dF/dx = 0, and 
F{x,fi) = Fo{iJ,). Fo(— /x) = Fq{ij,) because of electron reflection at the mirrors and the symmetry 
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of the loss cones. Up to first order, F{x,fi) = -Po(m) + and we have 



dFi _u d 
dx 2 dji 



(Al) 



We integrate this equation over x along a closed back and forth trajectory of a trapped electron 
shown by the dotted Unes in Figure 1(b), to obtain 



d/dfi [(1 - ^i'^)^Fo/^^l] + 2Fo/uTm = 0, 
Fo{-n) = Fo{n), Fo(±/Xcrit) = 0. 



(A2) 



We solve equation (A2) by a further expansion, l/vTm ^ 1- The even solution in the "inside" 
region 1 — |/x| S> e~^'^'^ up to first order is 



F^^ = C^^ 



1-J-ln ' 



, 1 - ImI » e" 



VTm 1 - 

On the other hand, the zero order solution in the "boundary" regions 1 — |//| <C 1 is 







C(^)ln 



1 — ^crit 



1 - yUcrit < 1 - < 1. 



(A3) 



(A4) 



We match solutions (A3) and (A4) together in regions e ^ 1 — <C 1 to finally obtain 
Tm = Inm, justifying l/uTm <^ 1. This is the first result in equation (5). 



B. Solution of equation (3) in the limits Agff <^ Im <^ A^/Agff and A^/Aeff <C Ir 



Let us consider the kinetic equation (3) in the more limited case X <^ 1^ (note that Aefr <C A) . 
This means that in the kinetic equation {V/iy){d/dx) < \/lm ^ 1- We will also show that 1/vTm ^ 
1. The validity of this assumption appears below. To zero order, we have dF/dji = 0, so F{x^ jj,) = 
Fq{x). Fq{—x) = Fq{x) because of symmetry. Up to first order, F{x,iJ,) = Fo{x) + Fi{x,iJ,), and 
we have 



2dn _ 



(1-M^ 



dFi 
dfx 



dx 







(Bl) 



We integrate the above equation over jj,, and then set = ±1 to find the constant of integra- 
tion. As a result, we obtain Fq/t„i <C V{dFQ/dx) [so, \/vTm is of second order], and dFi/dfi = 
— {V/i'){dFQ/dx). We integrate this last equation over fi once more, and obtain 



Fi = -{nV/i^)idFo/dx) + C{x), 
where C{x) is another integration constant. 



(B2) 



We continue the expansion of the kinetic equation (3) to next order. Up to second order, 
F{x,fi) = Fo{x) +Fi{x,^) + F2{x,fi). Using equation (B2), we have 



2dn 



(1-M^ 



dFo 



V ff^x dx 



Fh 



(B3) 
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We integrate equation (B3) over ^ from —1 to 1 and obtain 

d^Fo/d^x + {2,u/t^V^)Fo = 0, 



Fo{-x) = F^{x). ^^^^ 

Finally, we integrate this equation and obtain zero order solution for the time-dependent distribu- 
tion function (2) 

f{t, x) = e-^'^^Foix) = e-*/""- cos (x^/^^7J^^ ) , (B5) 

where we drop an unnecessary normalization constant of integration. 

Now, to find r^, we calculate the flux of escaping electrons through the two escape windows 
(see Figure 1) 

dN/dt = -2 j nVf{t, lm/2) dfi = -{V/m) e"*/"- cos [(Z„/2) VSi^/r^F^ j ^ (b6) 

/^crit 

where we use equation (B5) for f{t,lm/2). On the other hand, the flux is equal to the change of 
the total number of electrons 



/ {df/dt)dfidx = -{4/rm)e-'/^"^^TmVy3u sin (Z^/2) VSz^/r^F^ . (B7) 

-lm/2J-l L -I 



rlm/2 pi 

Equating the two formulas for dN/dt, we obtain 

(3/16)(i.r^/m2) = tan^ (^./si^l^J^^'^ . (B8) 

In the limit \ <^ Im ^ A^/Aeff the argument of the tangent above is small, so we expand 

the tangent and obtain Tm = Jy~^{lm/^cs)j while Fq ~ const. In the limit A^/Acff <^ Im the left 
hand side of equation (B8) is large, therefore, the argument of the tangent is 7r/2, and we have 
Tm = z^^^(3/'7r^)(/m/A)^ [the third line in equation (5)], i.e. the escape time is controlled by diffusion 
in x-space. In both limits l/vTm ^ 1, as we assumed above, (and of second order). 

Now, the limit Aeg ^ ^ A is still left. The result in this case is the same as the result in 
case X <^ Im <^ A^/Aeff. However, instead of solving the kinetic equation, we give the following 
qualitative arguments supported by our numerical simulations (see Figure 2). The relaxation time 
of the electron distribution in /x-space can be estimated as At^ ~ i/"^. The relaxation time in 
x-space can be estimated as the crossing time Atx ^ Im/V = ^~^{lm/^) in case Im ^ A, and as 
the time of diffusion across /Six ~ ^~^{lml^)^ in case A ^ Im- All relaxation times are small 
compared to the escape time t„i, i.e. At^, At^ ^ Tm for the entire range Aeff <^ Im ^ A^/Acff. This 
means that the distribution function is approximately constant in x and fi, say Fq 1, / e"*/^"". 
We then carry out calculations similar to those we used in formulas (B6) and (B7) to find that 
Tm = 2^~^(^m/Aeff) [the second line in equation (5)]. 
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C. The Additional electron flow produced by electric fleld 

Let us, for simplicity, consider the Lorentz gas. The results for the Spitzer gas are similar. 

First, we derive an estimate for the additional flow dJ- of electrons that are in an interval 
V € \y,V + dV) of the velocity space, produced by an electric field E due to the change of the two 
loss cones of a mirror trap. Let us consider only the principle mirrors, because they mainly control 
the diffusion of electrons (see Section 3). In this appendix, we denote their mirror strength (the 
principle mirror strength) as M . 

The principle mirror strength is of order of five, so in the case when the magnetic field decor- 
relation length is more than or approximately equal to the electron mean free path, Iq ^ \, the 
escape of electrons from the mirror trap is mainly controlled by their spatial diffusion, see Section 2. 
Thus, in this limit, the electrons "do not care" about the size of the loss cones, and therefore, no 
additional flow arises. 

In the case Iq ^ \ there is a non-zero additional flow dJ-. In Figure 1 (b) , because of the electric 
field, the loss cone on the left, //crit,-) is not equal to that on the right, /Xcrit,+ - The size of the two loss 
cones is estimated from the conservation of the electron magnetic moment, {1 — ijl^)V'^ / B = const, 
and from the conservation of energy, meF^/2 -|- eEx = const. We have 

/i^Ht,± ~ 1 - 1/M ± {eElu/meV^M), (CI) 

where Im is the spacing of the principle mirrors. 

Let dT^ and dT^ be the absolute values of the fluxes of the escaping electrons to the right 
and to the left respectively. Then, their sum is 

d^+ + dP^ = {Im/tm) 2nV^fodV, (C2) 

where 2'KV'^fodV is the number density of electrons expressed in terms of the Maxwellian zero 
order electron distribution function /q, and tm is the escape time, see equations (2), (6), (16) and 
Section 2. The actual electron flow, d^, is equal to the difference of and dJ^-, because they 
are in opposite directions. An estimate for the ratio dT /{dTj^ + dT-) is 



dT dJ^4- — d!F- 



dJ='+ + dT- dT+ + dT- 



l jidji — / ^djJL I I jidfj, + I 
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' Merit, + Merit,- 
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Now, using equations (C1)-(C3), we obtain the additional electron flow 

d^ « -27TieE/me)foilh/TM)dV. (C4) 

The factor {1\j/tm) in this equation is proportional to the spatial diffusivity of the electrons (pro- 
vided that the diffusivity is controlled by the principle mirrors, see Section 3). Thus, it is obviously 

1\j/tm = '^L^Li where Al and are the standard Lorentz mean free path and collision fre- 
quency, and i?D is the reduction of the spatial diffusivity reported in Section 3. Using that Al oc 1/^, 
fL oc V~^, and equations (16), (43), Vt = \/2kT/me, we finally obtain 

d^^-{l/2){2/'Kfl^ [k^/'^T^/'^Ej ml'^e^luK) Rnv^e-^'dv. (C5) 



Now we like to compare this result for the additional flow dJ- with the main flow dJ- produced 
by the electric field due to acceleration of particles. The later is 




m, 



v'^e-'^'dv, (C6) 



where we substituted function /i given by (17), and function S{v) = 'jESEiv) given by (33) and (39). 
As a result, 



Because the electrical current and the heat flow are mainly transported by superthermal electrons 
v'^ ~ 4, the additional flow produced by electric field due to non-equal loss cones can indeed be 
neglected in comparison with the main flow due to acceleration of electrons by electric field. 
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